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ABSTRACT 

We model the NFW potential to determine if, and under what conditions, the NFW halo appears 
consistent with the observed velocity fields of low surface brightness (LSB) galaxies. We present mock 
DensePak IFU velocity fields and rotation curves of axisymmetric and non-axisymmetric potentials 
that are well-matched to the spatial resolution and velocity range of our sample galaxies. We find 
that the DensePak IFU can accurately reconstruct the velocity field produced by an axisymmetric 
NFW potential and that a tilted-ring fitting program can successfully recover the corresponding NFW 
rotation curve. We also find that non-axisymmetric potentials with fixed axis ratios change only 
the normalization of the mock velocity fields and rotation curves and not their shape. The shape 
of the modeled NFW rotation curves does not reproduce the data: these potentials are unable to 
simultaneously bring the mock data at both small and large radii into agreement with observations. 
Indeed, to match the slow rise of LSB galaxy rotation curves, a specific viewing angle of the non- 
axisymmetric potential is required. For each of the simulated LSB galaxies, the observer's line-of-sight 
must be along the minor axis of the potential, an arrangement which is inconsistent with a random 
distribution of halo orientations on the sky. 

Subject headings: dark matter — galaxies: kinematics and dynamics 
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1. INTRODUCTION 

Cosmologically motivated numerical simulations of 
cold dark matter (CDM) describe very specifically the 
properties of the dark matter halos that should be ob- 
served in the universe. The simulations show that 
CDM halos are cuspy, meaning that the density of 
the halo, regardless of its mass, rises very steeply to- 



ward the center (e.g. Dubinski 1994 


;lNavarro et al.lll996l 


1997 


: Moore et al.lll999: Reed et al. 


20031 iDiemand et all 


2005 


). The simulations also dictate the range of permis- 



sible values of halo parameters based on the assumed 
cosmology of the simulations. The concentration c of a 
halo, for example, depends on the density of the universe 
at the time the halo forms, which in t urn depends on the 
adop ted values of h, fl m , erg, etc. ([Navarro et al.|[i~996l 
1997). Due to this intimate connection to cosmology, the 
values of halo parameters are not arbitrary. 

The most well-known description of CD M halo behav- 
ior is the cu spy NFW halo where p ~ r _1 (jNavarro et al.l 
119961 119971 ). The rotation curves of these halos are pa- 
rameterized by two numbers: the concentration, c, and 
a characteristic velocity, V2oo- These two parameters 
cannot freely vary, nor can they vary independently of 
the other; there is a correlation between c and VW) 
(e.g. iNavarro et aHll997t Uinell200d: iBullock et alll200lt 
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IWechsler et al.l [2002T ). This c — V200 relation, combined 
with the cosmological constraints on c, means that the 
expected rotation curve for a CDM halo of a given mass 
is well-determined. 

Though the need for dark matter in disk galaxies 
has long been in d icated by fla t rotation curves (e.g. 
iRubin et all 119801 : iBosmal Il98ll ). it has been less ob- 
vious that the dark matter halos are consistent with 
cuspy CDM halos. Because low surface brightness 
(LSB) galaxies are thought to be dark matter-dominated 
down to small radii (de Blok & McGaugh 1996; de 
Blok & McGaugh 1997; Bordello & Salucci 2001, but 
see Fuchs 2003), their kinematics have been used as 
probes of the density distribution of galaxy mass dark 
matter halos. Rotation curves derived from Hi veloc- 
ity fields and long-slit Ha observations are frequently 
consistent with halos having a cored p ~ r° density 
distribution (e.g. iFlores fe Primacklll994t Ide Blok et al.l 
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19961 Ide Blok fe Bosmal 120021: iMarchesini et al.l 120021; 
Cote. Carignan. fc Freemanl " I2000D rather than the 
steeper profile of the NFW halo. This result 
has also been supported by rotation curves de- 
rived from high-resolution two-dimensional velocity 
fields obtained with integral field spect r ographs (e.g. 
Chemin et ail 12004 iGentile et all 120051 : iSimon et al.l 
2005 HKuzio de Narav et al.ll2006l l2008D 

In iKuzio de Narav et alJ (|2006L 120081 hereafter K06 
and K08, respectively), we presented DensePak Integral 
Field Unit (IFU) Ha velocity fields, rotation curves, and 
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halo fits for a sample of 17 LSB galaxies. We fit both 
a cored pseudoisothermal halo (p ~ r°) and a cuspy 
NFW halo (p ~ r _1 ) to the data and found the halo 
central densities and rotation curve shapes to be better 
described by the cored halo model. The NFW fits to the 
DensePak rotation curves were often found to have con- 
centrations lower than what is expected for galaxies in a 
ACDM cosmology (see also Gentile et al. 2007, but see 
Swaters et al. 2003b for a different conclusion) and to fa- 
vor a power spectrum having a lower amplitude on small 
scales (Zentner & Bullock 2002, McGaugh et al. 2003, 
K08). We found that the NFW rotation curves speci- 
fied by the c — V200 relation (the rotation curves that 
our galaxies should have according to ACDM) are much 
more steeply rising than the observed DensePak rotation 
curves. In addition, these cosmologically consistent halos 
show a cusp mass excess at the centers of the galaxies, 
indicating that at least two times more mass is expected 
in the cuspy CDM halos than is allowed by the data. 

CDM halos must be both cuspy and follow the 
c — V200 rel ation defi ned b y ACDM. The density pro- 
files of the IK06I and IK08I data are not well-described 
by cuspy halos, nor do the galaxies fall on the c — 
V200 relation. These DensePak results are consistent 
with many p revious long-sl i t and Hi studies o f LSB 
galaxies (e.g. Ide Blok et all 1200 ll iBolatto etaTl 120021 : 
Ide Blok fc Bosmall2002l ISwaters et aT]l2003aD. as well as 
similar DensePak studies bv lSimon et alj (|2005t ). That 
different observational techniques (with different data re- 
duction and analysis procedures as well as sources of er- 
ror) lead to similar conclusions suggests that perhaps the 
discrepancy between the NFW halo and the observations 
does not arise at the telescope or during data analysis, 
but rather is due to an incorrect assumption about the 
specific form of the NFW halo potential. 

Our goal in this paper is to model the NFW halo to 
determine if, and under what conditions, it appears con- 
sistent with the observed data. Starting with an axisym- 
metric NFW potential, how must it be modified (e.g., 
introduction of an asymmetry) in order to appear con- 
sistent with both the observed two-dimensional velocity 
field and the derived rotation curve? We construct a 
model disk galaxy embedded in an NFW halo and then 
"observe" it in the same way as we have observed our 
sample of galaxies with DensePak. We then compare the 
mock velocity field and the derived mock rotation curve 
to the real galaxy data. 

We adopt a numerical approach to investigating 
non-axisymmetric halo potentials because once ax- 
isymmetry is broken, the data analysis becomes 
much more complicated. Noncircular motions and 
asymmetries are traditionally investigated by doing 
a higher-orde r Fourier decomposition of the veloc- 
ity field (e.g. Schoeiimakc rs. Franx. fc de Zeeuwl Il997t 
IWong. Blitz, fc Bosnia! 120041 ). We have tried this ap- 
proach, but it was not sufficiently well-constrained for 
these difficult LSB targets to give unique results. But 
we do find that useful constraints can still be extracted 
by simulating what is expected to be observed for various 
hypothesized halo potentials. 

The paper is organized as follows. In § 2 we describe 
the simulations. The axisymmetric NFW potential is ex- 
plored in § 3. In § 4 we describe the mock velocity fields 
and rotation curves produced by a non-axisymmetric 



TABLE 1 
Simulated NFW Halo Parameters 







V200 


^excess 


Galaxy 
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km s -1 


II 
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( 2 ) 


(3) 


(4) 


NGC 4395 


8.6 


87 


59 


DDO 64 


9.2 


62 


10 


UGC 4325 


6.9 


249 


40 


F583-1 


8.7 


83 


37 


F563-1 


8.4 


101 


13 


F583-4 


9.1 


67 


22 


UGC 5750 


9.1 


67 


31 


F563-V2 


7.9 


130 


22 


F568-3 


8.2 


110 


13 



Note. — Columns 2 and 3 list the NFW halo parameters for the 
simulated galaxies. Listed in column 4 are the radii at which the 
observed DensePak rotation curves and input NFW CO nstr rotation 
curves overlap in the minimum disk case. These are the radii out 
to which the observed and mock rotation curves are compared. 

NFW potential with a constant axis ratio. We deter- 
mine in § 5 the non-axisymmetric potentials that best 
describe the observed galaxy data. We discuss our re- 
sults and conclusions in § 6. 

2. DESCRIPTION OF SIMULATIONS 

TV-body simulations show that CDM hal os of all masses 
can be desc ribed by the NFW po tential (INavarro et al.l 
19961 119971) and i ts variants (e.g.. iDiemand et al I l2005t 



Moore et aij|1999t iNavarro et al.ll2004lReed et al.l 2003) . 
These cuspy halo potentials show a steep rise in the mass 
density toward the center of the halo. Most theoreti- 
cal estimates of the inner slope of the halo mass den- 
sity profile are as steep or steeper than that of NFW, 
so we choose to simulate the NFW halo as the conserva- 
tive case. If the NFW potential predicts a dark matter 
halo with a steeper density profile than is allowed by 
the observed galaxy data, then even more steeply rising 
potentials are automatically ex cluded. In some formula- 
tions (e.g.. INavarro et aill2004f ). there is no well-defined 
inner slope, which continues to roll over to a value that 
asymptotes to a flatter value than NFW. However, this is 
a small effect at small radii. Th e difference between t he 
original NFW profile and that of INavarro et alj (|2004D is 
too small to b e det ected observationally. 

In K06 and K08 we defined a constrained NFW halo, 
NFW cons t r . We required the halo to match the velocities 
at the outer radii of each galaxy by choosing a value of 
V200 which forced the NFW velocities to agree with the 
data points at large radii with the minimum requirement 
of falling within the e rrorbars of the data. We then used 
the c — V200 relation (jNavarro et al.l[T997t Ide Blok et al.l 
2003) to determine the corresponding cosmologically- 
consistent co ncentration. Th i s is ad justed to the 'vanilla' 
cosmology of iTegmark et alj (2004) by subtractin g 0.011 
dex in concentration (see [McGaugh "eFldl [2003h . Ac- 
cording to ACDM, these are the rotation curves that 
our galaxies should have. The chief remaining uncer- 
tainty in the normalization of the c — V200 relation is the 
power spectrum. Rot ation curves data prefer lower a% 
(|McGaugh et~aHl2057h . 

Our goal is to compare mock DensePak velocity fields 
and rotation curves of NFWconstr halos to the observed 
DensePak veloc ity fi elds and rotation curves of the LSB 
galaxies in IK06I and IK08I . We model those galaxies that 
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Fig. 1. — (a) Observed DensePak velocity field of UGC 4325. (b) Mock DensePak velocity field of the axisymmetric simulation. Both 
velocity fields have isovelocity contours at 10 km s _1 intervals, (c) Residual velocity field showing the differences between the UGC 4325 
data and the velocity field of the idealized (i.e. no velocity dispersion), axisymmetric NFW CO nstr halo, (d) Same as (c) but for the mock 
DensePak velocity field. The residuals are large and obvious in (c). Bottom: Observed and mock rotation curves. The solid points are 
the observed DensePak rotation curve of UGC 4325, the solid line is the NFW rotation curve corresponding to the input NFW potential, 
and the open (red) circles are the rotation curve recovered from the mock velocity field. The arrow indicates the radius out to which the 
rotation curves are compared. The last three points of the recovered mock rotation curve are high because of a lack of fibers at large radii. 



have well-sampled velocity fields and rotation curves that 
are constrained at large radii by previous long-slit and/or 
Hi rotation curves, allowing NFW cons t r halo fits to be 
made. Of our 17 galaxies, 9 meet these criteria. The 
spatial resolution and DensePak coverage of these data 
vary. The paramet ers o f the N FW eons t r halos in the limit 
of minimum disk (K06; IK08I ) for each of the 9 modeled 
galaxies are listed in Table 1. The galaxies are listed in 
order of decreasing spatial resolution, from NGC 4395 
(~ 20 pc/") to F568-3 (~ 375 pc/"). 

We developed a code which does a fourth-order Runge- 
Kutta (RK4) test particle integration of point masses 
moving in a two-dimensional rigid analytic NFW poten- 



tial. Specifically, the potential used is: 



GM 200 ln(l + -f ) 

Rf(c) ' ■ 



where 



and 



R = \A 2 + (y 2 /9 2 ), 

D R200 
c 

/(c) = ln(l+c) 



1 + c 



(1) 

(2) 
(3) 

(4) 



In these equations, M200 is the enclosed halo mass at ra- 
dius i?200j <7 is the axis ratio (q = y/x), and c is the con- 
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Fig. 2. — Same as Figure 1, but for the less well-resolved galaxy F583-4. 



centration of the halo (c = i?2oo/-Rs)- The halo param- 
eters are set to th ose of the NFWconstr halo determined 
for each galaxy in IK06I and IK08I and are listed in Table 
1. Each simulated galaxy is an infinitely thin exponen- 
tial disk of 10,000 test particles. For the nearby galaxies 
UGC 4325 and DDO 64, the number of test particles 
was increased to 100,000 to ensure sufficient sampling of 
particles to recreate the higher resolution data. The two- 
dimensional disk is in the plane of the potential. To wash 
out any numerical pattern noise of the initial conditions, 
we integrate for 50 half- mass rotation periods having 500 
timcstcps each. Each simulated galaxy is given the disk 
scale length, spatial resolution, and inclination of the real 
galaxy. 

Each simulated galaxy is then "observed" by 



DensePak. DensePak is an integral field spectrograph 
on the 3.5 m WIYN telescope at the Kitt Peak National 
Observatory (KPNO). It is a 43" x 28" fixed array of 3" 
fibers with 3.84" separations. We model the 85 work- 
ing fibers, as well as the 5 missing or broken fiber s, in 
the m ain bundle. For the galaxies observed in lK06l and 
IK08I . the fiber bundle orientation on the sky and the to- 
tal number of pointings per galaxy were tailored to each 
galaxy so that the critical central regions were covered 
by the DensePak fibers. We aim for obtaining roughly 
equivalent coverage of the simulated galaxies by using 
similar numbers and alignments of DensePak pointings 
on the simulations. Prior to extracting a rotation curve, 
these mock DensePak velocity fields are given the veloc- 
ity dispersion observed in the real DensePak galaxy ve- 
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locity fields. We have defined the velocity dispersion of 
the DensePak data to be the fiber-to-fiber velocity varia- 
tion; to recreate this in the mock data, we randomly add 
the desired dispersion to the fibers in the mock velocity 
field. Rotation curves were then derived from t he mock 
Dens ePak observations by using the NEM O (jTeubenl 
fl995l) program ROTCUR (|Begemanl fl989l) . ROTCUR 
treats the observed velocity field as an ensemble of tilted 
rings and fits for the center, systemic velocity, inclina- 
tion, position angle, an d rota tion velo city in each ring. 
The reader is referred to IK06I and IK08I for a more exten- 
sive explanation of ROTCUR and its application to the 
DensePak velocity fields. 

3. AXISYMMETRIC NFW HALOS 

The most obvious and simple starting point is to as- 
sume an axisymmetric halo potential. The axis ratio q is 
equal to 1 and the test particles move on circular orbits. 
With this straightforward potential, we can test whether 
or not DensePak observations are sufficient to detect the 
signature of NFW halos in the velocity fields and/or 
whether the data analysis procedure with ROTCUR also 
suffices to recover NFW rotation curves. 

In Figure 1 we model the NFW constr halo of UGC 4325 
and "observe" the simulation with 5 pointings of the 
DensePak array. The pointings are arranged to match 
the spatial coverage of the real galaxy as much as pos- 
sible. Observed and residual velocity fields for both 
UGC 4325 and the simulation are shown. The simulated 
galaxy has the same fiber-to-fiber velocity dispersion as 
the real galaxy: a = 9.0 km s _1 . UGC 4325 is one of the 
most nearby (D rj 10 Mpc) and well-resolved galaxies in 
our sample. Diffuse Ha emission was abundant in the 
galaxy and was detected in almost all of the DensePak 
fibers. 

Figure 1 demonstrates two important points. First, 
as evidenced by the very small residuals of the mock 
DensePak velocity field, the DensePak instrument is able 
to successfully detect an NFW velocity field. The resid- 
uals are generally < 5 km s _1 across the entire observed 
area. This means that observed DensePak velocity fields 
are not inconsistent with NFW halos because of an inad- 
equacy of the experimental design or analysis. Second, 
the observed DensePak velocity field of UGC 4325 is not 
consistent with the axisymmetric NFW constr halo; most 
of the residuals are ~ 10 km s^ 1 , and there is a signifi- 
cant region of ~ 15 km s" 1 residuals near the center. For 
UGC 4325, and the other LSB galaxies in our sample, 15 
km s _1 residuals are non-trivial. The observ ed fib er-to- 
fiber velocity dispersions are ~ 6-10 km s _1 (KOff), and 
since mass scales as a 2 , the implied mass difference is a 
factor of two or more. In addition, noncircular motions 
caused by disk instabilities, such as spiral or bar modes, 
are expected to be small in LSB galaxies. The low sur- 
face mass densities of the disks provide little self-gravity 
to drive such modes, and their high dark matter con- 
tent provides a higher degr ee of stabilization than in high 
surfa ce brightness galaxies (|Mihos. McGaugh. fc de Blold 
[1993). 

At the bottom of Figure 1 are the observed and mock 
rotation curves derived from the velocity fields. The re- 
covered mock rotation curve is consistent with the in- 
put rotation curve (x 2 . = 0.93) out to R exC ess ~ 40", 
where R excess is defined to be the radius at which the 




Observer's Line-of-Sight 



Fig. 3. — Orientation (rf>) of the elongated axis of the two- 
dimensional non-axisymmetric potential with respect to the ob- 
server's line-of-sight. In the <j> = 0° case, the elongated axis of 
the potential is along the observer's line-of-sight, whereas in the 
<p = 90° case, it is perpendicular to the observer's line-of-sight. 
[See the electronic edition of the Journal for a color version of 
this figure.] 

observed DensePak rotation curve and the input axisym- 
metric NFW cons t r rotation curve begin to overlap. This 
shows that accurate extraction of the rotation curve of an 
axisymmetric NFW halo with ROTCUR is also possible. 

In Figure 2 we show similar plots for F583-4. This 
galaxy has lower spatial resolution (D w 49 Mpc) than 
UGC 4325 and only a single pointing of DensePak cov- 
erage. The velocity fields, residuals, and rotation curves 
show that despite the reduced sampling and lower data 
quality, an axisymmetric NFW halo can still be detected 
if present. While the differences between the observed 
and mock velocity field residuals are not as pronounced 
as in the UGC 4325 case, the input NFW rotation curve 
is successfully recovered by ROTCUR (jfi = 0.32). 

The velocity field and rotation curve data and simu- 
lations plotted in Figures 1 and 2 together show that 
the DensePak IFU and the tilted-ring fitting program 
ROTCUR are able to successfully identify an axisymmet- 
ric NFW halo in dat a of b oth high and low quality, if one 
is present. That the lK06l and lK08l samples of DensePak 
observations are inconsistent with NFW halos suggests 
that if the underlying halo potential is NFW, it must 
not be an axisymmetric NFW potential. This is perhaps 
not surprising, as CDM si mulations suggest that the halo 
poten tials are triaxial (e.g. lHavashi. Navarro. &; Springe! 
12001 . 

4. NON-AXISYMMETRIC NFW HALOS WITH A FIXED 
AXIS RATIO 

We next consider non-axisymmetric two-dimensional 
NFW potentials with axis ratios q < 1 that are con- 
stant with radius. These 2D potentials are equivalent to 
3D prolate dark matter halos in which the long axis of 
the halo coincides with the elongated axis of the disk. 
We simulate halos with axis ratios q = 0.98, 0.96, 0.94, 
0.92, 0.90, 0.88, 0.86, and 0.84, similar to the range 
of non-axisymmetry seen in the CDM simulations of 
lHayashi et al.l (|2007f ). Because axisymmetry has been 
broken (q ^ 1), the test paticles are no longer moving 
on circular orbits and not all lines of sight in the plane 
of the disk are equivalent. This means that the observed 
mock DensePak velocity field and derived rotation curve 
are affected not only by the value of q, but also by the 
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Fig. 4. — Mock DensePak velocity fields for the non-axisymmetric NFW simulations of UGC 4325. UGC 4325 is a well-resolved galaxy 
with multiple pointings of DensePak coverage. Simulations with an axis ratio q = 0.98 are shown in the left column, and simulations with 
q = 0.88 are in the right column. The observer's veiwing angle changes from = 0° in the top panels to <j> = 90° in the bottom panels. 
Note that the effect of viewing angle is much more pronounced in the q = 0.88 case. For easy comparison, all of the velocity fields are on 
the same color /velocity scale, and isovelocity contours are drawn at 10 km s — 1 intervals. 



orientation (0) of the potential's elongation with respect 
to the observer's line-of-sight (see Figure 3). The poten- 
tial is elongated along the observer's line-of-sight in the 
= 0° case, whereas in the cf> = 90° case, the potential 
is elongated perpendicular to the observer's line-of-sight. 
For 0° < cj) < 90°, the elongation is at an intermediate 
viewing orientation. For each value of q, the orientation 
of the potential is set to cf> = 0°, 30°, 45°, 60°, and 90°. 

4.1. "Observed" Mock DensePak Velocity Fields 

In Figures 4 and 5 we show a series of mock DensePak 
velocity fields for two galaxies representative of the range 
of data presented in |K0l and lK08t UGC 4325 and 
UGC 5750. As previously mentioned, UGC 4325 has 
high spatial resolution and extended DensePak coverage. 
As in the axisymmetric case, 5 pointings of the DensePak 
array are overlayed on these new non-axisymmetric sim- 
ulations. In contrast, UGC 5750 is a more distant galaxy 



(D w 56 Mpc) and both the real galaxy and the simula- 
tions have only one pointing of DensePak coverage. 

These are simulations of non-axisymmetric NFW ha- 
los that obey the cosmic c — V200 relation. The virial 
velocity V200 has be en ch osen to m atch each galaxy (the 
NFW co „ str . halos of IK06I and IK08D . In this section, we 
explore the effect of introducing a non-axisymmetric po- 
tential with equal squashing q at all radii. 

In the q = 0.98 simulations, the potential is nearly cir- 
cular. Throughout their orbits, the particles maintain a 
roughly constant distance from the center of the poten- 
tial and as a result, have approximately constant orbital 
speeds. The viewing angle therefore has little effect on 
the observed velocity field. The q = 0.98, = 0° -> 90° 
mock velocity fields appear very similar, looking not only 
to be consistent with different realizations of the same 
underlying potential, but also very much like the mock 
velocity field of the axisymmetric potential. 
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Fig. 5. — Same as Figure 4, but for the mock DensePak velocity fields for the non-axisymmetric NFW simulations of UGC 5750. 
UGC 5750 is a more distant galaxy and has only a single pointing of DensePak coverage. 



The same cannot be said for the mock velocity fields of 
the q — 0.88 simulations. With orbits deviating signifi- 
cantly from circular, a particle's orbital speed depends on 
its location, making the viewing angle quite important. 
In the = 0° orientation, the particles moving along the 
observer's line-of-sight are traveling along the long axis 
of the potential and are moving at the maximum orbital 
speed. These particles are moving faster than particles 
on circular orbits at the same radius. The minimum- 
maximum velocity range observed by DensePak is larger 
than what is observed in the axisymmetric case, and the 
derived rotation curves will reach higher velocities. The 
opposite situation is happening in the (ft = 90° orienta- 
tion. In this case, the particles moving along the ob- 
server's line-of-sight are traveling at the minimum or- 
bital speed and the minimum-maximum velocity range 
observed by DensePak is smaller than what is observed 
in the axisymmetric case. The rotation curves derived 
from these data will therefore be suppressed. The differ- 
ence between the observed velocity ranges of the (ft = 0° 
and (ft = 90° velocity fields becomes more exaggerated 



the more noncircular the potential becomes. 

Some (q, (ft) combinations can automatically be ex- 
cluded as possible descriptions of the observed DensePak 
galaxy data based simply on the mock velocity fields they 
produce. The observed DensePak galaxy velocity fields 
put constraints on the allowable velocity range of the 
mock velocity fields, as well as the correlation between 
velocity and position. Regardless of how the mock ro- 
tation curves may turn out, if the observed and mock 
velocity fields do not match, the corresponding simula- 
tion is not a viable solution. For example, the (g, (ft) — 
(0.88, 0°) mock velocity field of UGC 4325 shown in Fig- 
ure 4 can rule out that particular axis ratio/viewing ori- 
entation combination for that galaxy. Overall, the mock 
velocity field covers a much larger velocity range than 
the UGC 4325 data, and when the velocities at the same 
positions in the two velocity fields are compared, they are 
inconsistent over a large portion of the observed area. 

Because of the rapidly rising velocities at the centers 
of NFW halos, the isovelocity contours of NFW velocity 
fields are pinched in the central regions (jde Blok et alJ 
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Fig. 6. — A comparison of three mock velocity fields to the ob- 
served UGC 4325 data using velocities measured along two slits 
placed parallel to, and offset from, the minor axis of the veloc- 
ity fields. The solid (green) line is for the (q, <f>) = (0.88, 0°) 
mock velocity field, the long-dash (red) line is for the axisym- 
metric mock velocity field, the short-dash (blue) line is for the 
(q, </>) = (0.88, 90°) mock velocity field, and the line+circles are 
the UGC 4325 data. At the same distance from the center of the 
velocity field, the (q,<j>) = (0.88, 90°) mock velocity field has ve- 
locities most similar to the observed data. A typical errorbar is 
shown in the lower left corner. [See the electronic edition of the 
Journal for a color version of this figure.] 

I2003D . This pinch is a distinctive signature of the cuspy 
NFW halo. If the velocity field is noisy or has high ve- 
locity dispersion, the pinch is more difficult to see. We 
can quantify the pinch by measuring velocities along slits 
that arc offset from, and parallel to, the minor axis of the 
velocity field. It is in this fashion that we compare the 
appearance of the UGC 4325 DensePak velocity field and 
three mock velocity fields. In Figure 6, we have placed 
a 3" x 28" slit (the width of a DensePak fiber and the 
width of the DensePak array, respectively) 6" away from 
each side of the minor axis (about the separation of two 
rows of DensePak fibers) of the observed DensePak veloc- 
ity field of UGC 4325 and the axisymmetric NFW mock 
velocity field (both shown in Figure 1), as well as the 
(q, 4>) = (0.88, 0° and 90°) mock velocity fields in Figure 
4. When we plot the average of the measured velocities 
as a function of position along the slit, we can readily 
see that of the 3 mock velocity fields, it is the one pro- 
duced by the (q, </>) = (0.88, 90°) potential that is most 
like the data. At the same position in the velocity fields, 
the axisymmetric, and especially the (q, <f) — (0.88, 0°), 
mock observations are detecting velocities much higher 
than the galaxy data. 

In the next section we derive mock rotation curves for 
all the mock DensePak velocity fields and examine the 
effect that the asymmetry of the potential has had on 
both the normalization and shape of the derived rotation 
curves. 

4.2. Derived Mock Rotation Curves 

Th e LSB galaxies observed with DensePak in lK06l and 
IK08I surely have some level of noncircular motions, but 
we have assumed only circular motion when deriving 
the rotation curves with ROTCUR; we treat the mock 
observations the same way. The test particles are, by 



construction, no longer on circular orbits in these non- 
axisymmetric potentials, but because we did not correct 
for this in the real data, we do not correct for it in the 
mock observations. Any errors in the galaxy rotation 
curves which may have resulted from the assumption of 
circular motion will be reproduced in the rotation curves 
of the mock observations. 

In Figure 7 we show the mock rotation curves for the 
non-axisymmetric simulations of the well-resolved galaxy 
UGC 4325. Each panel shows the observed rotation 
curve of UGC 4325 along with the <j) = 0°, 30°, 45°, 
60°, and 90° mock rotation curves for a single value of q. 
Similar to the trends seen in the mock velocity fields, we 
find that as the axis ratio decreases and the potential be- 
comes increasingly more elongated, the influence of the 
viewing angle on the inferred rotation curves becomes 
more striking: the mock rotation curves for the different 
values of spread farther out in velocity space. We can 
quantitatively measure and compare the shapes of these 
rotation curves using the ratio of the radii containing 80% 
and 50% of the velocity at R eX cess- The average values 
of R 8 o/R 5 o for the UGC 4325 q = 0.98, 0.90, and 0.86 
mock rotation curves, for example, are 2.7±0.7, 2.8±0.8, 
and 3.1±0.7, respectively. This indicates that the over- 
all shapes of the mock rotation curves are not changing 
significantly as q or <f> change; it is the normalization of 
V(r), including V max , that is shifting up or down. This is 
an important point to recognize, as it means that simply 
adopting a different (lower) value of V200 f° r the under- 
lying NFW halo (which cannot be done without ignoring 
cosmological constraints) will not reconcile the observed 
and mock rotation curves; the mock NFW rotation curve 
will still not fit the data properly. Though none of the 
mock rotation curves match the shape of the entire ob- 
served rotation curve of UGC 4325, those that are most 
consistent with the data at small radii are the <f> = 90° 
mock rotation curves. This is true for all of the values 
of q that were simulated, but it is the (f> = 90° rotation 
curves in the q < 0.90 simulations that have the most 
overlap with the data. 

As a comparison to the well-resolved observations of 
UGC 4325, we show in Figure 8 similar plots of the mock 
rotation curves for UGC 5750. Despite the lower spatial 
resolution, we find the mock rotation curves of UGC 5750 
to behave in very much the same way as the mock ro- 
tation curves of UGC 4325. We again see that the am- 
plitude, not the shape, of the rotation curve changes as 
q and (f> change, with the magnitude of the change be- 
coming more pronounced as the potential becomes more 
asymmetric. We also find that the cj> = 90° mock ro- 
tation curves are again the closest to approaching the 
data, though substantial overlap at small and intermedi- 
ate radii does not occur until q < 0.86. 

We examine in Figure 9 how changing the input NFW 
halo parameters affects the mock rotation curves, specif- 
ically exploring if q and <fr can change the shape of a 
more slowly rising input NFW rotation curve. For both 
UGC 4325 and UGC 5750, we have simulated NFW ha- 
los with low values of V200 and the lowest corresponding 
concentration within the scatter of iBullock et alj (|2001 ) 
(UGC 4325: c = 5.2, V 2Q0 = 140; UGC 5750: c = 6.9, 
V200 = 60). These new input NFW rotation curves fall 
between the rotation curves representing the -la scatter 
on c and V200 of the original NFW CO nstr rotation curves. 
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Fig. 7. — Mock rotation curves (lines) for the non-axisymmetric NFW simulations of the nearby galaxy UGC 4325 (points). Each panel 
is for a different value of the axis ratio q (labeled in the lower right). The solid lines are for = 0°, 30°, 45°, and 60°. The dashed line is 
for <f> = 90°. For all values of q, the <f> = 90° line has the most overlap with the data at small radii. The dotted line in the lower right panel 
is the input NFW CO nstr rotation curve. The arrows indicate the radius out to which the rotation curves are compared. [See the electronic 
edition of the Journal for a color version of this figure.] 



In addition, these slowly rising rotation curves overlap 
the observed DensePak data at small radii, in contrast 
to the NFW conS (r rotation curves which match the data 
at large radii. Essentially, we choose to match the inner 
rather than outer velocities with halos drawn from the 
favorable edge of the plausible cosmological distribution. 

For both galaxies, we find the new mock rotation 
curves to behave similarly to the mock rotation curves 
in Figures 7 and 8. As the axis ratio q decreases, the 
mock rotation curves for the different values of <j> scatter 
about the input NFW rotation curve, spreading farther 
out in velocity space. Even though these new mock rota- 
tion curves overlap some of the observed DensePak data 
at small radii, none are formally acceptable fits (Xr ^ 
1). More importantly, the shapes of the mock rotation 
curves are not significantly changing as q and 4> change: 
the average values of Rso /R50 for the UGC 4325 q = 0.98, 
0.90, and 0.84 mock rotation curves are 2.3±0.3, 2.8±0.4, 



and 3.1±1.2. The mock UGC 5750 rotation curves are 
essentially flat between ~10" and ~30", preventing use- 
ful measurements of Rso /Rso- From Figure 9 we can see 
that regardless of the input NFW halo parameters, q and 
4> change only the normalization, not the radial behavior, 
of the mock rotation curves. 

It is also worth stressing that, given the behavior of the 
mock rotation curves in Figures 7 and 8, as well as Fig- 
ure 9, observers should see rotation curves with a range 
of normalizations: there should be rotation curves both 
above and below the nominal rotation curve expected 
from the c — V200 relation (compare the mock rotation 
curves to the dotted lines in Figure 9 and the lower right 
panels of Figures 7 and 8). But this is, in fact, not what is 
observed in long-slit data (iMcGaugh. Rubin. &: de Blold 
[200lt IdeBlok fc Bosmal 12002ft. LSB and NFW rotation 
curves nearly always differ in the sense that the observed 
rotation curve velocities at small radii must increase so 
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Fig. 8. — Same as Figure 7, but for UGC 5750. [See the electronic edition of the Journal for a color version of this figure.] 



that the data match the models, or equivalently, the 
NFW rotation curve velocities must decrease so the mod- 
els match the data. LSB rotation curves which are higher 
than NFW rotation curves are seldom found, if ever. 

There is a trade-off between q and (f> such that dif- 
ferent combinations of the two parameters can produce 
similar mock rotation curves. In the following section, we 
explore what combination of (q, </>) minimizes the differ- 
ences between the NFW halo and the DensePak galaxy 
observations. We simulate the NFW conS ( r halos rather 
than low V200 halos like those in Figure 9 because the 
constrained halos were required to match the velocities 
at the outer radii of each galaxy, a reasonable constraint 
since dark matter must explain the high velocities at 
large radii where the contribution of the baryons has 
fallen off. Parameter space is too large to explore all 
possible initial halos. However, given that plausible com- 
binations of NFW c and V200 parameters give rather de- 
generate rotation curves, and that q and qb affect only the 
normalization, not the shape of V(R), our choice should 
lead to fairly general results. 



5. MINIMIZING THE CUSP MASS EXCESS WITH Q AND 4, 

In EH we showed that there is a substantial cusp 
mass excess near the centers of the galaxies when the 
NFWconstr halo is used to describe the dark matter halo. 
Evaluating the difference between the NFW CO nstr rota- 
tion curve and the observed galaxy rotation curve in 
terms of mass rather than velocity, we determined that 
interior to the radius where the two rotation curves be- 
gin to overlap {Rexcess), NFW halos are at least twice 
as massive as the galaxy data will allow. In this sec- 
tion, we are interested in determining for each DensePak 
galaxy what combination of (q, (f>) minimizes the differ- 
ences between the observed and mock rotation curves 
out to R excess where, in the limit of zero stellar mass, 
the cusp mass excess is ~ (see Table 1). In this fash- 
ion, one can imagine a toy model in which the halo of a 
particular galaxy is squashed to the best fit (q, <p) within 
Rexcess while outside of Rexcess we have a more nearly 
spherical, cosmologically consistent NFW co „ S ( r halo. 

Figures 7 and 8 showed that regardless of the axis ra- 
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Fig. 9. — Top Ro w: The solid line is the NFW CO nstr rotation curve. The short dashed lines are the NFW rotation curves corresponding 
to the ±1<t scatter {Bullock ct al. 2001) expected in the concentration in ACDM, and the long dashed lines correspond to the ±lcr scatter 
on V200- The halo of this galaxy could plausibly be drawn from anywhere in this range. For example, the dotted line is the NFW rotation 
curve representing a low V200 and the lowest corresponding concentration within the scatter. This halo provides a good initial match to 
the inner data at the expense of falling well short of the outer data. Lower panels: Mock rotation curves for the non-axisymmetric NFW 
simulations of this low V200 halo (dotted lines). Each panel is for a different value of the axis ratio q (labeled in the upper left). The 
solid lines are for = 0°, 30°, 45°, 60°, and 90°. The arrows indicate the radius out to which the rotation curves are compared. [See the 
electronic edition of the Journal for a color version of this figure] 



tio, the mock 4> — 90° rotation curves came closest to the 
observed rotation curves of UGC 4325 and UGC 5750. 
To confirm that the differences between the rotation 
curves derived from the simulations and from the ob- 
served galaxy data are really minimized at <f> ps 90° and 
not somewhere between qb = 60° and 90°, we ran addi- 
tional simulations at <j> = 75°, 85°, 86°, 87°, 88°, and 
89°. We then determined for each combination of q and 
(f> how well, as measured by x% the mock and observed 
rotation curves matched out to R eX cess- 

In Figure 10, we plot the best <f> for each value of q for 
each galaxy. For nearly all of the 9 simulated galaxies, 
the mock rotation curves are the most consistent with 
the DensePak galaxy rotation curves when <j» is between 
85° and 90° for all values of q. A <j) — > 90° means that the 
elongated axis of the NFW potential is pointing perpen- 



dicular to our line-of-sight. This required 4> is completely 
inconsistent with a random distribution of halo orienta- 
tions on the sky. 

There is one galaxy in our sample, F563-V2, which has 
a preferre d value of (j) other than 90°. This galaxy has a 
bar in it (jPildis et al.lll997h . We ran additional simula- 
tions for F563-V2 found that <j> 1=3 55° is the optimal view- 
ing angle for matching the mock rotation curves to the 
observed rotation curves. This position angle matches 
that of the bar. It would therefore seem that we have de- 
tected the expected noncircular motion associated with 
the bar rather than the sq uashing of the halo (see also 
ISpekkens fc Sellwoodl[2007f) . 

This result for F563-V2 confirms that we are able to 
detect the presence and orientation of an asymmetry in a 
velocity field. If the other DensePak galaxies contain bars 
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Fig. 10. — The viewing angle cj> that minimizes the differences be- 
tween the mock and observed galaxy rotation curves out to Rexcess 
for each value of q for each galaxy. Nearly all of the galaxies fall 
on top of each other in the <j> = 85° to 90° range for all values of q. 
The optimal viewing angle for F563-V2 (open squares) is <f> ~ 55°; 
there is a stellar bar at this position angle. The results for F563-1 
and NGC 4395 are not shown because <j> is unconstrained due to 
the radial extent of the data. Nevertheless, the fact that we detect 
the bar in F563-V2 is an encouraging confirmation of the method. 
[See the electronic edition of the Journal for a color version of this 
figure] 
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Fig. 11. — The value of Xr as a function of <f> for the best-fitting 
q for each galaxy. Although the <j> — » 90° mock rotation curves are 
not formally good fits to the observed data, they are the minimum 
in x 2 -space. The exception is F563-V2, for which Xr ~ 1 f° r the 
optimal viewing angle of <f> « 55° . [See the electronic edition of the 
Journal for a color version of this figure.] 

or are embedded in non-spherical NFW halos, we would 
be able to detect the asymmetry. That tj> rs 90° for all 
the other DensePak galaxies demonstrates that either the 
data are inconsistent with non-axisymmetric NFW halos, 
or that we must accept the unlikely coincidence that all 
of these galaxies are oriented such that the elongated axis 
of the potential is perpendicular to our line-of-sight. It 
is not surprising that the effect goes in this sense as the 
rotation curves of LSB galaxies are persistently measured 



to be shallower than expected for NFW halos. 

There are two galaxies whose results are not shown 
in Figure 10: F563-1 and NGC 4395. The value of 
is unconstrained for both of these galaxies due to the 
radial extent of the data. For F563-1, there are only 
a few data points to compare between the observed and 
mock rotation curves. NGC 4395 is a very nearby galaxy 
(D « 3.5 Mpc), and although there are many data points 
to compare in the rotation curves, the data probe a radius 
of less than ^800 pc. 

As discussed in § 4.2, q and 4> can be used to change 
the amplitude of the NFW rotation curve. They do not, 
however, alter the overall shape of that rotation curve. 
This is reflected by high Xr values for the comparisons 
of the mock and observed rotation curves. Although the 
Xr values are typically greater than 1, the sharp decline 
in x r as 4> 90° indicates that <p ~ 90° is truly the 
minimum (see Figure 11), even though the mock rotation 
curves are not formally acceptable fits to the observed 
data. In Figure 12 we plot the "best-fitting" -> 90° 
rotation curves (cb — > 55° for F563-V2) over the observed 
galaxy data for each galaxy. The two galaxies with the 
highest spatial resolution, DDO 64 and UGC 4325, are 
clear examples of how the mock rotation curve has shifted 
down in velocity such that the inner half of the mock 
rotation curve is roughly consistent with the observed 
data, but the outer half of the mock rotation curve falls 
below the observed data. 

Despite not being able to fully match the observed 
galaxy rotation curves within R exceS s, the important 
trend in Figures 10 and 11, and reinforced by Figures 4-6, 
is that (j> is being driven toward 90° if one wants to match 
the observed data at small radii where the cusp/core 
problem is most severe. This means that the elongated 
axis of the NFW potentials for every DensePak galaxy 
(with the exception of the barred galaxy F563-V2) must 
point perpendicular to our line-of-sight. This required 
<j) is completely inconsistent with a random distribution 
of halo orientations on the sky. A non-axisymmetric po- 
tential with a fixed axis ratio may be able to bring parts 
of the NFW rotation curve into agreement with the ob- 
served data for individual galaxies, but in general, very 
peculiar, observer-dependent conditions must occur. 

It is worth mentioning here that the mismatches be- 
tween the "best-fitting" mock rotation curves and the 
data, as seen in Figure 12 for example, are a result of 
real differences in the velocity fields; information is not 
being lost or suppressed as ROTCUR collapses all the 
data contained in the two-dimensional velocity fields into 
a one-dimensional representation of the rotation. As was 
shown in Figures 4 and 5, the lowest velocity portions 
of a particle's orbit are being preferentially observed in 
the <j) — 90° viewing orientation, ensuring that the ob- 
served range of velocities detected in the mock velocity 
fields is both small and slowly varying. But like the ro- 
tation curves show, only some parts and not all of the 
observed and mock data can be made to match. Shown 
as an example in Figure 13 are the observed UGC 4325 
DensePak velocity field, the "best-fitting" <j> = 90° ve- 
locity field, and the residual velocity field showing the 
differences between the two. While the residuals are rel- 
atively small in the central regions of the velocity field, 
there are multiple areas at larger radii where the resid- 
uals are in the range of ~ 10-15 km s _1 , or more. As 
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Fig. 12. — The "best-fitting" <j> — > 90° rotation curves (tf> — > 55° for F563-V2) for each galaxy. The solid (red) line is the average of the 
<j> = 85° — » 90° (<j> = 45° — * 60°) rotation curves and the shaded (blue) band outlines the spread in those rotation curves. For comparison, 
the dotted line is one realization of the cf> = 60° (0 = 90°) rotation curve. The arrow indicates Rexcess- From top to bottom, the galaxies 
are ordered by increasing distance. The simulations of UGC 4325 do not extend to the outermost observed rotation curve point. The inner 
7" of the mock rotation curves of UGC 5750 are poorly sampled, as in the real data. The results for F563-1 and NGC 4395 are not shown 
because <f> is unconstrained. [See the electronic edition of the Journal for a color version of this figure.] 



was discussed in § 3, in galaxies where the fiber-to-fiber 
velocity dispersions are measured to be ~ 6-10 km s , 
15km s" 1 residuals imply a significant mismatch. 

6. DISCUSSION AND CONCLUSIONS 

In this paper we have simulated the two-dimensional 
NFW halo and tested several modifications to the po- 
tential in an attempt to simultaneously reconcile both 
the NFW velocity field and rotation curve with observed 
DensePak galaxy data. We have found that it is dif- 
ficult to make the cuspy NFW halo appear consistent 
with core-like data without violating the predicted range 
of NFW parameters expected in the ACDM cosmological 
model. 



Beginning with simulations of an axisymmetric NFW 
potential, we found that both the DensePak IFU instru- 
ment and the rotation curve fitting program ROTCUR 
are able to successfully identify the NFW potential. Ob- 
served galaxy data is inconsistent with the velocity fields 
and rotation curves corresponding to axisymmetric NFW 
potentials. The simulated observations show that our 
data would detect the NFW cusp if it were present. 

We also tried a non-axisymmetric potential with a fixed 
axis ratio. We found that if the parameters of the NFW 
halo (determined from the c — V200 relation) are held 
constant and only the axis ratio and viewing orientation 
are varied, parts of the mock velocity fields and mock 
rotation curves, but not their entire area or length, could 
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Fig. 13.— (a) Observed DensePak velocity field of UGC 4325. 
(b) The "best-fitting" non-axisymmetric <j> = 90° mock NFW ve- 
locity field of UGC 4325. (c) Residual velocity field showing the 
difference between the observed and mock velocity fields. The ve- 
locity fields are similar in the central regions, but become more 
mismatched at large radii. 

be made to roughly match the observed galaxy data. The 
axis ratio and viewing orientation work to change only 
the normalization, not the radial behavior, of the mock 
data. The shape of the predicted NFW rotation curve 
remains distinct from the observations. This remains 
true even if a more slowly rising NFW rotation curve is 
simulated. 

Only when the elongated axis of the mock galaxies 
is oriented perpendicular to the observer's line-of-sight 
(4> ~ 90°) are the critical velocities at small radii in both 
the velocity fields and rotation curves consistent with 
observations. Having all halos elongated perpendicular 
to our line of sight is clearly not reasonable. This con- 
straint on (f> can be relaxed if we start with an NFW 
halo having a lower V200 (see Figure 9), but changing 
the halo parameters requires that we disregard cosmolog- 
ical constraints. Numerical simulations have shown that 
scatter in the c — V 200 relation exists (A(log c) = 0.18; 
iBullock et "all ()2001h ). but this is not enough to "fix" the 
mock NFW data, as the values of the halo parameters 
would have to be outside the range of the allowable dis- 
persion. Furthermore, even though lowering V200 allows 
for more scatter in the problem with the shape of the 
mock NFW rotation curve remains unresolved. 

In order to reconcile both the entire area of the NFW 
mock velocity fields and the entire length of their derived 
rotation curves with galaxy data, we need an asymme- 
try that preferentially suppresses velocities at small radii. 
The asymmetric NFW potentials with fixed axis ratios 
that we have tested in this paper have altered the ve- 
locities at all radii, either suppressing all the measured 
velocities (the 4> — 90° case) or boosting all of them (the 
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<f) —* 0° cases). One possible way to address that prob- 
lem ma y be to invoke a non- a xisymmetry tha t varie s with 
radius (jHavashi et al.ll2007h . lHavashi etall (|2007[ ) have 
suggested that galaxy-sized CDM halos are triaxial with 
radially varying axis ratios. They find the halo poten- 
tial to be highly elongated near the center (b/a — > 0.78 
and c/a — > 0.72) and increasingly more spherical at large 
radii. This is the general behavior that the results of our 
two-dimensional simulations suggest is required. 

It remains to be seen though if such an asymmetry is a 
viable solution to the problem. The simulations will al- 
ways be constrained by the fact that LSB galaxy velocity 
fields and rotation curves are slowly rising. This will lead 
to the problem of a preferential viewing angle that we 
have already encountered. Any potential that deviates 
significantly from axisymmetry at the center will have 
to be viewed at an angle which lowers, not increases, the 
observed velocities; the inner ellipsoid will need to be per- 
pendicular to the line-of-sight. On the other hand, if one 
insists that a radially-varying asymmetry is the correct 
solution and that the viewing orientation is randomly 
distributed, both slowly and rapidly rising LSB galaxy 
rotation curves should be observed. LSB galaxy rota- 
tion curves that are steeper than NFW rotation curves 
are not generally found. Of t he 50+ long-slit rota- 
tion curves in the literature (e.g . Zackrisson et al. 120061; 
Spekkens. Giovanelli. fe Havnesll2005t Ide Blok fc Bosmal 
2002tlMcGaugh. Rubin, fe de Blokll2001h . it is very com- 
mon to see slowly rising rotation curves, but exceedingly 
rare to see rotation curves that are in excess of the ex- 
pectation for NFW halos that obey the ACDM c - V200 
relation. 

The analysis of an NFW potential with a variable axis 
ratio is sufficiently complex that if one is to do it right, 
it should be done in three dimensions with interacting 
particles. If having a highly elongated potential is key 
to "fixing" the NFW velocity field and rotation curve, 
then adiabatic contraction should also b e considered, as 
it will round out t he inner halo potential (|Dubinskilll994h 
Bailin et al. 2007). Adiabatic contraction also has the ef- 
fect of increasing the dark matter density of the halo over 
its initial val ue and serves to worsen the concentration 
probl em (e.g.. lGnedin et al.l[200llSellwood fe McGaughl 
|2005| ) . We therefore defer the analysis of a radially vary- 
ing axis ratio to a forthcoming paper in which the com- 
plexity of the simulations is increased by moving to three 
dimensions and including gas physics. 
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